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SECTION  I 


INTRODUCTION 


1.  Sub  ject  Matter 

The  determination  of  the  aerodynamic  characteristics  of  missiles 
and  projectiles  from  free-flight  ballistic  range  data  is  dependent  on 
three  main  factors:  (1)  accurate  ballistic  range  orientation  and  position 
data,  (2)  appropriate  analytical  aerodynamic  models  for  use  in  parameter 
estimation  programs  which  adequately  account  for  configuration  and  flow 
field  characteristics,  and  (3)  estimation  algorithms  which  provide 
reliable  estimates  of  the  aerodynamic  coefficients  using  the  given  data 
and  model.  This  report  deals  with  a study  of  methods  to  improve  the 
accuracy  of  the  second  of  these  important  factors. 

The  emphasis  of  the  study  is  placed  on  the  development  of  a model 
description  procedure  which  properly  correlates  flight  item  configuration 
and  kinematics  with  major  flow  field  characteristics.  The  motivation 
behind  this  particular  approach  is  to  provide  for  the  selection  of  a 
model  form  on  the  basis  of  configuration  and  expected  flow  field  parame- 
ters. 

2.  Historical  Background 

Since  the  least  squares  extraction  method  of  Chapman  and  Kirk 
(Reference  1)  came  into  use  in  the  analysis  of  ballistic  range  data, 
increasingly  complex  forms  of  aerodynamic  models  have  been  employed  in 
various  parameter  estimation  programs.  All  of  these  models  have  been 
based  on  the  series  expansion  of  aerodynamic  forces  and  moments  in  terms 
of  various  kinematic  variables.  The  coefficients  in  these  expansions, 
usually  called  stability  derivatives,  are  then  determined  by  the  parameter 
estimation  process.  The  recent  enlargement  of  the  viable  envelope  for 
ballistic  testing  brought  about  by  improvements  in  experimental  and 
analytical  methods  has  put  considerable  emphasis  on  the  need  for  improved 
aerodynamic  models  as  well. 

Two  indirect  methods  for  model  structuring  have  been  recently 
reported.  In  Reference  2 the  Estimation  Before  Modeling  (EBM)  method 
is  described.  This  technique  has  the  particular  advantage  that  a state- 
dependent  model  need  not  be  selected  prior  to  state  identification.  Thus, 
total  force  and  moment  values  are  determined,  and  particular  model  forms 
for  comparison  with  data  from  other  sources  can  be  determined  over  very 
restricted  regions  of  the  parameter  space  (e.g.,  over  a specified  range 
of  angles  of  attack).  Reference  3 describes  recent  advances  in  model 
identification  methods.  The  general  technique  here  is  to  analytically 
describe  an  aerodynamic  model  which  is  very  general  in  terms  of  independent 
variables  which  can  be  used  to  describe  the  forces  and  moments.  Then 


the  particular  model  is  selected  based  on  comparisons  with  range  data. 

This  method  thus  determines  the  model  form  and  the  coefficient  values 
through  Interaction  between  the  model  form  and  the  extraction  algorithm. 
The  method  allows  for  change  in  model  form  when  the  physical  behavior 
of  the  system  within  a given  range  of  independent  variables  requires  that 
such  a change  be  made.  In  this  way  the  method  indirectly  Incorporates 
a correlation  with  the  physical  flow  field  behavior. 

Preliminary  studies  of  aerodynamic  models  were  performed  by  the 
Principal  Investigator  during  a USAF/ASEE  Summer  Research  Fellowship 
at  Eglin  AFB.  These  studies,  which  concentrated  on  determining  exist- 
ing analysis  methods  and  Initiating  mass-offset  studies  and  investiga- 
tions into  model  form  are  presented  in  Reference  4.  The  results  of 
that  Investigation  were  sufficiently  encouraging  to  lead  to  a more 
detailed  formulation  of  the  model  selection  process.  The  present  report 
gives  the  details  of  the  investigation  procedure  which  was  developed. 

3.  Scope 

The  remainder  of  this  report  is  devoted  to  an  explanation  of  the 
studies  which  led  to  the  formulation  of  a research  plan  for  the  improve- 
ment of  aerodynamic  models  for  use  in  ballistic  range  data  analysis. 

These  studies  consist  of  three  major  phases:  (1)  review  of  current 
analytical  methods  in  the  field  of  parameter  estimation  to  determine 
model-dependency,  particularly  with  regard  to  analytical  form, 

(2)  establishment  of  a data  base  from  which  aerodynamic  model  inferences 
may  be  made  and  (3)  analytical  determination  of  mass-offset  effects  on 
dynamic  motion  through  derivation  of  various  transformation  relations. 

The  development  of  the  analytical  tools  required  to  implement  these  model 
studies  is  also  carried  out. 

Section  II  contains  a description  of  the  analytical  schemes  which 
are  currently  being  used  by  various  agencies  to  determine  the  aerodynamics 
of  flight  items  through  analysis  of  free-flight  data.  The  emphasis  of 
this  discussion  is  concerned  with  the  role  of  the  aerodynamic  model  in 
the  extraction  process. 

Section  III  Identifies  the  type  and  source  of  aerodynamic  data 
for  three  general  classes  of  flight  item.  These  are:  the  axisymmetrlc 
body  alone,  the  axisymmetrlc  body  with  tail  surfaces  and  the  axisymmetrlc 
body  with  wing  and  tall  surfaces.  Special  attention  is  paid  to  non-linear 
aerodynamic  characteristics,  since  the  current  difficulties  with  model 
selection  are  predominantly  because  of  these  non-linear  terms.  Also 
included  in  this  section  is  a discussion  of  the  adaptability  of  several 
different  model  forms  in  the  coefficient  extraction  process. 

Section  IV  provides  the  details  of  the  general  coordinate,  inertia 
and  aerodynamic  coefficient  transformations  necessary  to  effect  a mass- 
offset  analysis  of  a free-flight  test  item.  Sensitivity  equations  are 
also  derived  which  indicate  the  degree  to  which  a basic  parameter  is 
affected  by  an  axis  offset. 
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Section  V outlines  a specific  research  plan  incorporating  the 
results  of  this  study.  The  purpose  of  this  plan  is  to  provide  a 
systematic  technique  for  evaluating  a variety  of  aerodynamic  model 
forms.  Block  diagrams  of  the  model  evaluation  procedure  are  presented 
together  with  proposed  techniques  for  providing  a correlation  between 
aerodynamic  model  formulation  and  identifiable  flow  field  characteris- 
tics. 


Section  VI  is  a summary  of  the  information  in  this  report.  Important 
conclusions  and  special  features  of  the  study  are  presented  here. 

In  Appendix  A the  essential  features  of  the  mass-offset  transformations 
are  presented.  In  Appendix  B program  listings  for  simulation  and 
transformation  calculations  are  presented. 
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PARAMETER  ESTIMATION  METHODS 
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1 . Introduction 

There  are  four  methods  of  parameter  estimation  which  have  direct 
application  to  the  free-flight  data  analysis  problem  and  which  appear 
to  be  at  the  forefront  according  to  recent  references.  The  details  of 
these  methods  are  important  in  the  present  work  only  with  regard  to  the 
role  which  the  aerodynamic  model  plays  in  the  computations.  Particular 
attention  is  paid  to  the  effect  of  various  model  forms  on  the  estimation 
calculations , 

The  basic  methods  under  consideration  here  are;  the  extended 
Kalman  filter  (EKF),  maximum  likelihood  method  (MLM) , estimation  before 
modeling  (EBM)  and  model  structure  determination  (MSD) . One  other  method 
was  given  careful  consideration  during  the  evaluation  period.  The  i 

square-root  variable  metric  technique  described  in  Reference  5 was  ruled  ' 

out  as  a possible  method  for  ballistic  analysis  because  the  numerical  | 

derivative  computations  require  more  data  for  accurate  implementation 
than  is  usually  available  from  ballistic  tests. 

There  is  some  evidence  to  suggest  that  the  form  of  the  aerodynamic 
model  affects  the  accuracy  of  the  extraction  computations  if  an  unduly 
large  number  of  coefficients  are  being  extracted.  This  tendency  precludes 
representation  of  the  aerodynamic  terms  in  very  general  terms.  It  is  | 

important  to  be  able  to  Identify  important  physical  effects  which  are 
expected  in  the  flight  envelope  and  to  reflect  these  in  the  model  so 
that  the  number  of  cefflclents  being  sought  can  be  as  small  as  possible. 

Based  on  evidence  from  numerous  applications  of  parameter  estimation 
methods  to  ballistic  data  analysis  at  Eglin  AFB,  the  order  in  which 

coefficients  are  extracted  also  affects  the  accuracy  of  the  outcome.  ' 

This  factor  was  considered  in  the  study  of  the  extraction  methods  | 

discussed  here. 

[ 

Finally,  it  should  be  emphasized  that  the  purpose  of  this  section 
is  not  to  select  an  extraction  algorithm,  but  rather  to  determine  any  I 

special  features  which  a model  should  have  to  be  effective  in  a particular  r 

method.  Currently,  the  maximum  likelihood  method  is  being  used  for  ' t*, 

data  analysis  at  the  Aeroballlstic  Research  Facility  at  Eglin  AFB.  Further  [• 

model  studies  at  that  installation  will  be  made  with  that  extraction  1 

algorithm. 

2.  Extended  Kalman  Filter  P 

The  Kalman  filter  theory  determines  the  state  of  a system  by 
utilizing  measurements  of  some  or  all  of  the  system  state  variables  i 

to  make  estimates  of  the  system  condition  at  the  next  interval. 

Actually,  the  original  formulation  applied  only  to  linear  systems, 
but  the  EKF  method  can  be  applied  to  highly  non-linear  systems.  This 
method  of  analysis  is  a special  case  of  the  more  general  maximum 
likelihood  method. 
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The  application  of  this  method  to  the  free-flight  trajectory 
analysis  problem  follows  the  same  fundamental  steps  as  the  original 
Chapman-Kirk  method.  That  is,  an  analytical  trajectory  is  computed 
using  a specified  aerodynamic  model.  This  analytical  model  is  then 
compared  with  discrete  trajectory  parameters  measured  in  the  ballistic 
range.  A least  squares  minimization  is  then  performed  yielding  a set 
of  aerodynamic  coefficients  which  provides  the  best  analytical  fit  to 
each  experimental  data  point  in  the  trajectory.  The  basic  difference 
between  the  two  methods  lies  in  the  minimization  algorithm.  The  EKF 
algorithm  permits  solutions  of  the  coupled  dynamic  equations  while 
the  Chapman-Kirk  method  requires  that  the  translational  equations  be 
solved  Independently  from  the  rotation  equations.  Both  methods, 
however,  commonly  employ  polynomial  expansions  to  represent  the 
aerodynamic  terms  in  the  equations  of  motion.  This  model  format  is 
common  to  most  dynamic  performance,  stability  and  control  studies 
(see  Section  V),  and  it  simplifies  correlation  with  experimental 
results  from  other  forms  of  t-  ?t. 

There  is  no  inherent  restriction  in  the  EKF  method  on  the  model 
form  as  long  as  the  analytical  formulation  can  be  put  in  terms  of  the 
state  variables.  Some  difficulties  have  been  reported,  however,  in  the 
operation  of  the  algorithm  if  initialization  is  not  carefully  achieved. 

The  EKF  performance  is  very  sensitive  to  the  Initialization  of  parameters, 
so  a model  which  cannot  be  well  evaluated  initially  may  lead  to  a 
poor  filter  performance  even  though  the  model  itself  may  be  very  good. 

3.  Maximum  Likelihood 

This  method  is  a more  general  modification  of  the  basic  least 
squares  method  of  Chapman  and  Kirk  which  utilizes  a different  optimi- 
zation scheme  to  achieve  coefficient  extractions.  In  that  method  the 
solutions  were  achieved  for  the  uncoupled  equations  of  motion 
because  the  scales  of  the  deviation  matrices  for  translational  and 
rotational  motions  were  different,  precluding  the  coupled  analysis. 

The  maximum  likelihood  algorithm  applies  weighting  factors  to  the 
traslational  and  rotational  deviation  functions  so  that  the  best  fit  of 
the  analytical  model  to  the  experimental  data  can  be  achieved  by  maxi- 
mizing a likelihood  function  (that  the  estimate  and  the  actual  data 
coincide) . The  aerodynamic  model  plays  the  same  role  in  this  model  as 
it  does  in  the  Chapman-Kirk  and  EKF  methods.  There  seems  to  be  nothing 
Inherent  in  the  method  which  would  rule  out  any  model  form  with 
undetermined  coefficients  for  describing  the  aerodynamic  forces  and 
moments.  However,  it  is  essential  that  the  formulation  be  posed  in  such 
a way  that  the  equations  of  motion  can  be  analytically  differentiated 
with  respect  to  the  coefficients  so  that  the  sensitivity  equations  used 
in  the  statistical  analysis  can  be  determined. 

This  technique  is  currently  used  at  Eglin  AFB  to  determine 
aerodynamic  coefficients  from  ballistic  data.  One  disadvantage 
in  doing  model  studies  with  this  or  the  method  previously  discussed 
is  that  analytical  expressions  are  required  for  the  evaluation  of 
sensitivity  parameters.  Thus  an  Internal  change  in  the  computer  code 
must  be  made  with  each  model  change.  The  form  of  the  model  could,  thus, 
affect  its  Implementation  into  the  extraction  program  because  of 
coding  or  storage  difficulties,  but  the  quality  of  the  estimation 
should  not  be  compromised  because  of  this  difficulty. 
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Estimation  Before  Modellni 


This  technique  is  quite  different  from  those  previously 
described  in  that  an  aerodynamic  model  is  not  required  prior  to 
the  first  estimation  calculation.  Furthermore,  the  results  of 
this  analysis  are  in  the  time,  rather  that  state,  domain,  so  the 
total  aerodynamic  forces  and  moments  are  determined  only  as  time 
functions. 

It  is  assumed  that  the  motion  of  the  flight  item  is  described 
by  a system  of  non-linear  differential  equations  (the  equations  of 
motion)  in  which  the  state-dependence  of  certain  terms  (e.g., 
aerodynamic  forces  and  moments)  is  unknown.  These  unknown  terms  are 
replaced  by  special  spline  time  series  with  unknown  coefficients. 

These  coefficients  are  then  determined  by  conventional  means  (e.g., 
maximum  likelihood)  through  comparison  with  experimental  data.  Once 
the  time  domain  solutions  become  available,  the  total  aerodynamics 
are  known.  That  is,  the  forces  along  each  axis,  the  moments  about 
each  axis  and  any  auxiliary  state-dependent  functions,  such  as 
control  forces  and  moments,  become  completely  determined  in  the  time 
domain.  Correlation  of  these  results  with  the  state-time  relationship 
provided  by  the  equations  of  motion  yields  the  values  of  the  total 
aerodynamic  characteristics  for  any  given  value  of  a fundamental 
state  variable  or  for  a value  of  a state-varlable-determlned  quantity 
(such  as  angle  of  attack  or  Mach  number) . 

This  technique  provides  the  total  aerodynamic  force  and  moment 
as  a function  of  time.  However,  no  detailed  information  regarding 
the  contributing  parts  of  the  total  is  provided.  Thus,  the  resulting 
total  provides  little  direct  information  regarding  the  effects  of 
variation  of  parameters  on  the  total  result.  This  makes  the  result 
impossible  to  check  directly  with  other  test  data  (e.g.,  wind  tunnel 
data).  In  addition,  the  total  force  and  moment  functions  are  totally 
dependent  on  the  degree  of  fit  achieved  by  the  estimation  process. 

One  facet  of  this  application  which  is  appealing  as  a possible 
modeling  technique  is  the  use  of  splines  for  function  approximation 
in  the  analysis.  The  details  of  various  spline  functions,  as  outlined 
om  Reference  6,  Indicate  that  a complex  function  can  be  described  with 
only  a few  spline  regions  provided  the  approximated  functions  are 
smooth.  It  is  felt  that  this  approach  is  worth  additional  study  with 
regard  to  applications  to  conventional  stability  derivative  functions 
which  seem  not  to  be  adequately  described  by  standard  polynomial  functions, 

5,  Model  Structure  Determination 

A very  recent  innovation  in  the  broad  field  of  parameter  esti- 
mation and  system  performance  evaluation  is  the  application  of  statis- 
tical data  comparisons  to  the  determination  of  the  basic  structure  of 
the  models  which  Influence  system  behavior.  As  applied  to  ballistic 
range  data  analysis  this  method,  called  model  structure  determination 
(HSD) , would  repetitively  attempt  to  match  experimental  data  with  a 
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variety  of  models  for  the  aerodynamic  forces  and  moments.  The  basic 
extraction  algorithms  could  be  Identical  to  those  described  In  parts 
2 and  3 earlier,  but  the  analysis  and  statistical  testing  goes  beyond  that. 
By  this  technique  any  number  of  model  forms,  each  with  Its  own  set  of 
undetermined  coefficients,  can  be  Included  In  the  search  for  a model 
or  group  of  models  which  will  provide  the  best  simulation  of  the 
experimental  trajectory.  Such  an  approach  allows  for  the  simulation 
of  rapidly  fluctuating  aerodynamic  terms  which  may  arise  from  such 
things  as  vortex  roll-up,  shock  stall  or  wake  Impingement. 

The  general  model  form  required  In  this  analysis  Is  one  for 
which  appropriate  sensitivity  relations  can  be  derived  by  differenti- 
ation with  respect  to  the  parameters.  This  obviously  makes  the  choice 
very  wide,  ranging  from  a simple  algebraic  polynomial  to  a Fourier 
series.  The  number  of  different  model  forms  and  the  generality  of 
each  Is  limited  by  the  computational  constraints  of  time  and  storage. 

In  any  given  case  the  model  forms  are  selected  based  on  the  physical 
phenomena  they  are  expected  to  describe. 

6.  Summary 

The  studies  which  have  been  discussed  here  have  been  based 
completely  on  Information  available  In  the  open  literature  and  on 
actual  computational  results  which  have  been  achieved  by  the  Ballistics 
Branch  at  Eglin  AFB,  Florida.  The  Principal  Investigator  has  only 
used  the  Chapman-Klrk  and  maximum  likelihood  methods  with  polynomial 
aerodynamic  models  as  they  are  formulated  in  the  operational  programs 
at  Eglin  AFB.  Nevertheless,  some  conclusions  of  a general  nature  can 
be  concluded  from  the  present  study  of  model/extraction  method  Interaction. 

For  a great  many  flight  situations  the  aerodynamics  can  be 
accurately  represented  by  conventional  algebraic  polynomial  expansions. 

In  all  of  the  methods  described  the  polynomial  can  be  conveniently 
Incorporated.  Furthermore,  the  "stability  derivative"  format  is 
easily  compared  with  results  from  other  types  of  tests.  However,, 
there  are  some  physical  phenomena  which  cannot  be  simulated  with  this 
model.  These  Include  forces  and  moments  In  regimes  where  neither 
even  nor  odd  function  characteristics  are  evident  because  of  complex 
Interaction,  In  flow  fields  containing  significant  vortex  growth  and 
for  cases  In  which  violent  maneuvers  are  encountered  because  of  surface 
stall.  Other  series  type  formulations  could  be  used  to  simulate  these 
cases.  For  example,  transcendental  series,  splines  or  any  of  the 
special  form  series  such  as  Fourier  or  Legendre  could  be  used. 

It  appears  Imperative  in  the  continuation  of  these  model  studies 
that  unusual  physical  phenomena  which  are  observed  to  occur  In  the 
flight  dynamics  of  various  configurations  must  be  correlated  to  specific 
model  types  so  that  a systematic  model  development  program  can  evolve. 

Of  course,  such  a study  requires  experimental  data  on  a variety  of 
configurations  and  flight  conditions.  The  next  section  describes 
Initial  Investigations  In  this  category. 


SECTION  III 


AERODYNAMIC  DATA  COLLECTION 


1.  Introduction 


In  order  to  develop  a systematic  approach  to  the  model 
formulation  task  it  is  necessary  to  establish  some  development 
guidelines.  One  of  these  guidelines  Involves  configuration  sel- 
ection. The  motivation  behind  establishing  certain  classes  of 
configuration  types  is  to  simplify  the  task  of  identifying  those 
parts  of  the  aerodynamic  model  which  are  strongly  configuration- 
dependent. 

The  present  research  plan  is  predicated  on  defining  models 
for  three  classes  of  flight  item.  These  are:  bodies  of  revolution 
without  control  surfaces,  bodies  of  revolution  with  tail  fins  only 
and  bodies  of  revolution  with  tall  and  wing  surfaces.  Of  course,  a 
large  volume  of  wind  tunnel  data  are  available  in  the  technical  lit- 
erature of  the  past  thirty  years  on  configurations  of  all  types, 
including  the  three  classes  of  Interest  here.  The  emphasis  in  seek- 
ing data,  however,  has  been  placed  on  establishing  modeling  criteria 
in  the  non-linear  regions  of  the  flight  profile  such  as  at  high  angle 
of  incidence  or  in  regions  where  Induced  effects  are  important.  Typi- 
cal data  sources  and  types  for  each  of  the  three  configuration  classes 
will  now  be  given. 

2.  Body  Alone 

Probably  more  data  collection  has  been  done  on  bodies  of 
revolution  than  on  any  other  class  of  flight  item.  The  data  ranges 
from  basic,  drag  and  moment  data  to  induced  forces  and  moments  for 
spinning  bodies.  The  computations  to  determine  the  linearized  character- 
istics (low  angle  of  incidence)  at  all  speeds  are  outlined  in  such  basic 
sources  as  References  7,8  and  9.  In  addition  NACA  and  NASA  have  pub- 
lished numerous  reports  of  wind  tunnel  tests  of  a variety  of  body  confi- 
gurations at  all  flight  speeds.  As  long  as  the  data  are  at  moderate 
angle  of  attack  they  can  be  adequately  modeled  by  conventional 
polynomials.  Indeed  such  characteristics  as  lift,  drag  and  moment 
variations  with  angle  of  attack  and  Mach  number  hold  no  surprises. 

Magnus  forces  and  moments  and  other  incidence  and  rotation  in- 
duced effects  are  not  so  well  known.  Reference  10  contains  data 
on  cone-cylinder  and  ogive-cylinder  which  would  be  difficult  to  fit  with 
a regular  polynomial  at  high  angles  of  attack.  However,  the  data  do 
appear  to  be  amenable  to  polynomial  fitting  on  a piecewise  basis. 
Reference  11  contains  fundamental  pressure  data  for  six  nose-cyllnder- 
boattall  bodies  of  revolution  for  angles  of  attack  to  60<’  at  super- 
sonic Mach  numbers.  Integration  of  these  data  should  yield  sufficient 
Information  to  model  some  of  the  high  angle  non-linearities  evident 
in  force  and  moment  data.  In  fact,  when  used  in  conjunction  with  the 
distributed  vortex  model  of  Reference  12,  a picture  of  vortex  effect 
modeling  should  evolve. 
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A very  important  part  of  the  data  base  is  a large  collection  of 
ballistic  data  accumulated  by  the  AFATL  staff  at  Eglln  AFB,  A 
variety  of  projectile  configurations  have  been  tested  balllstlcally 
and  in  the  wind  tunnel  to  provide  a very  good  set  of  related  data.  Most 
of  these  configurations  are  of  the  projectile  type,  however,  rather 
than  being  of  a shape  usually  associated  with  a missile  or  airplane. 

Also,  they  are  necessarily  spin  stabilized,  so  the  aerodynamic  behavior 
at  high  angle  of  attack  cannot  simulate  that  for  a non-spinning  body. 

In  light  of  the  experience  of  many  years  of  projectile  testing,  it  is  no 
surprise  that  the  ballistic  data  can  be  fit  with  very  good  accuracy 
by  the  aerodynamic  models  currently  in  use  (Reference  13).  In  fact, 
a recent  investigation  (Reference  lA)  indicates  that  accurate  estimates 
can  be  achieved  even  when  the  projectile  contains  a moving  Internal 
component . 

3.  Body-Tall  Fin 

The  most  common  type  of  guided  missile  is  the  finned  body. 

There  are  many  configurations  of  finned  missiles  which  consist  of  a 
body  of  revolution  with  tail  surfaces.  These  can  be  spin  or  fin  stab- 
ilized. The  fin  stabilized  missile  Introduces  two  classes  of  aerodynamic 
non-linearity  that  are  not  usually  encountered  in  the  unflnned  geo- 
metries. These  are  fin-body  Interference  (and  vise  versa)  and  body 
vortex  shedding  at  Incidence.  In  addition,  the  fin  geometry,  position 
on  the  body  and  body  geometry  all  Influence  the  aerodynamics.  Fortunate- 
ly there  has  been  a considerable  amount  of  testing  associated  with  finned 
missiles.  A configuration  known  as  the  Basic  Flnner,  which  consists 
of  a cone-cylinder  body  with  four  rectangular  planform  tail  fins  in 
horizontal  and  vertical  planes,  has  been  tested  extensively  in  free 
flight  and  in  the  wind  tunnel  (e.g.,  see  Reference  15).  Supporting 
analytical  work  for  configurations  at  low  angles  of  attack  is  also 
available  from  a variety  of  sources  such  as  Reference  16. 


Availability  of  data  for  these  configurations  at  high  angle  of 
incidence  is  less  secure.  A recent  report  of  a ballistic  research 
project  at  Eglln  AFB  is  encouraging  (Reference  17).  Extensive  wind 
tunnel  and  free  flight  tests  have  been  performed  on  an  ogive-cylinder 
body  with  swept  leading  edge  tall  fins.  Angles  of  attack  on  some  of  the 
tests  are  near  30°  and  distance  non-linearities  are  observed.  It 
appears  that  the  aerodynamic,  models  could  be  extended  with  a matched 
polynomial  (like  a spline)  to  Incorporate  these  areas  of  interest. 

Visual  studies  of  the  flow  have  not  been  made  to  determine  the  probable 
cause  for  the  high-angle  behavior.  Similar  studies  have  been  per- 
formed by  NASA  (Reference  18)  on  a finned  body  with  much  lower  fine- 
ness ratio  and  swept  vane-type  fins. 

A computer  search  of  government  documents  has  failed  to  turn 
up  test  data  for  tall-finned  configurations  at  high  incidence  published 
since  1970.  Most  of  the  basic  testing  on  these  configurations  was 
probably  performed  in  the  1950' s and  is  available  in  the  government 
archives.  On  the  other  hand,  the  fin  stabilized  missiles  are  not 
intended  to  operate  in  the  high  incidence  regime  and  extensive  testing 
in  that  regime  may  never  have  been  authorized. 
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The  model  forms  which  are  currently  being  employed  for 
studying  conventional  configurations  at  low  angle  of  incidence 
appear  to  be  adequate.  These  involve  representation  of  the  forces 
and  moments  as  polynomial  functions  of  Mach  number,  sine  of  the  total 
incidence  angle  and  the  aerodynamic  roll  angle. 

4 . Body-Tail-Wing 

The  advent  of  new  high  performance,  high  maneuverability 
aircraft  during  this  decade  has  necessitated  the  development  of  new 
offensive  and  defensive  missile  weaponry.  Missiles  now  in  the  inventory 
frequently  fly  in  a hlgh-g,  high  incidence  angle  environment.  These 
configurations  represent  the  ultimate  challenge  to  ballistic  testing 
and  data  analysis.  It  is  with  these  configurations  that  the 
greatest  innovations  in  model  formulation  will  be  necessary.  Typical 
of  the  data  with  severe  non-linearities  which  result  from  such 
configurations  are  the  test  results  shown  in  Reference  19.  This 
particular  configuration  experiences  severe  non-linear  excursions  in 
weathervane  and  lateral  stability  at  angle  of  attack.  Proper  modeling 
of  just  this  one  effect  would  require  a radical  departure  from  the 
conventional  polynomial  representation.  Similar  non-linear  behavior 
in  the  longitudinal  mode  is  reported  in  Reference  20  for  recent 
maneuvering  missile  configurations.  It  may  well  be  impossible  to 
isolate  each  physical  effect  which  is  causing  the  extreme  non-linearities 
because  the  configurations  are  complex  and  the  flow  fields  very  difficult 
to  analyze. 
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SECTION  IV 


MASS-OFFSET  CALCULATIONS 


1.  Introduction 

In  ballistic  testing  it  is  important  to  know  the  exact  center  of 
gravity  (c.g.)  location  of  the  model  because  the  free-fllght  behavior 
will  be  strongly  influenced  by  this  location.  Model  measurements  prior 
to  testing  determine  the  longitudinal  position  of  the  c.g.  and  the 
moments  of  inertia  of  the  body  about  geometrical  axes  of  symmetry 
through  this  longitudinal  c.g.  position.  The  data  reduction  of  the 
test  flight  is  then  predicated  on  the  c.g.  lying  on  the  axis  of  symmetry. 
If  the  c.g.  is  displaced  from  this  axis,  it  is  necessary  to  analyze 
the  dynamics  based  on  that  displaced  position.  Such  an  analysis  requires 
that  the  corresponding  inertias  and  aerodynamic  coefficients  also  be 
transformed  to  the  new  position  so  that  the  Newtonian  equations  of 
motion  will  apply. 

The  fundamental  concepts  of  these  transformations  were  given  in 
a very  general  form  in  Reference  4.  A summary  of  this  work  plus  some 
newer  developments  on  coefficients  and  sensitivities  is  given  in  the 
following  paragraphs. 

2.  Coordinate  Transformations 


Figure  1 illustrates  the  basic  coordinate  relationships.  Using 
the  angular  and  translational  displacements  as  parameters,  the  coordinates 
are  transformed  by 

' LBf  [X^l 

The  matrix  elements  for  these  and  the  following  transformations  are 
shown  in  Appendix  A. 

3.  Inertia  Transformations 


The  Inertia  tranformations  between  c.g.  axes  and  axes  of  symmetry 
are  given  by 


= [C]  [ij  + <•''-2*) 
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4. 


Velocity  Transformations 


The  translational  and  rotational  velocity  components  tranform 
as  follows: 


[UJ  - 

[B]  fu.-Vi 

(IV-3a) 

lu,]  = 

(IV-3b) 

II 

(IV-4a) 

- 

(IV-4b) 

5.  Force  and  Moment  Transformations 

The  total  forces  and  moments  can  be  transformed  according  to  the 


following  relations: 

II 

(IV-5a) 

II 

LU 

[QVih] 

(IV-5b) 

= 

[BUm.-TT)] 

(IV- 6a) 

ll 

(IV-6b) 

6.  Aerodynamic  Coefficient  Transformations 

In  order  to  apply  the  Newtonian  equations  of  motion,  the  axis 
system  must  be  located  at  the  c.g.,  and  the  forces  and  moments  must  be 
defined  with  respect  to  these  axes.  In  the  usual  case  the  c.g.  is  on 
the  axis  of  symmetry,  and  the  simulation  can  proceed  normally.  In  the 
case  of  a flight  item  with  mass-offset,  however,  the  axes  are  displaced, 
and  corrections  must  be  made.  These  corrections  Include  aerodynamic 
coefficient  transformations. 

For  most  flight  items  the  axis  displacements  will  be  very  small, 
and  approximate  corrections  can  be  made.  However,  in  this  report  a 
general  formulation  is  made  without  any  small  displacement  approxima- 
tions. A typical  derivation  of  one  of  the  aerodynamic  coefficients  is 
shown  in  the  following  paragraph.  Other  coefficients  can  be  derived 
in  a similar  manner.  A listing  of  the  key  coefficients  is  shown  in 
Appendix  A, 
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Ci)nsider  the  case  where  characteristics  known  with  respect  to  a^es 
{x2,  yg,  Zp);  an  aerodynamic  reference  ajtes  set,  are  to  he  transformed  to 
the  center  of  gravity  axes  set  (x,  , y^  , ).  Typically, 


- (bn  Fx2  t b?l  -f  bsi  ) 5./, 

But  typically  4 4- 

Fran  (3)  and  (^4)  these  derivatives  can  he  evaluated  to  he 

L "^Pc' 

DuT,  D33  j 'Du/,  ■ D22  ’j  'Sur,  “ O 


Thus  1^'  = 

Note  also  that 

Now,  let 


• sincxz  = 

VojCoSo/, 

9cXz  ' 

"j  ^ cos  0^2 

3 P>^z 

■DtX^  uTj.  ” '^0 

III 

0 

[(5?V«0A  ] 

(IV-7)‘ 

(IV-7a) 

(IV-7b) 

(IV-7c) 

(IV-7d) 

(IV-7e) 


Suhstitutian  of  equations  (7a)  through  (7e)  hack  into  equation  (7)  yields, 
for  a case  where  certain  derivatives  are  zero  (e.g.  dFy^/BuT,  =0 

~ [c^S^cos€cos(r  Cx^  4-  2 (co^.<^sin?6co5V4  S\o<Fcose  Sin2(r)G  ] 

' ^ *^2 

-/^2[2C(55^Sin26  CXu  + (coS<5sm^£a)5(r+  Sin<rSineSin(r)C2  , ] (IV-7f) 

'-*2.  U9 


where 


\ _ . 

Ai  C0S0I2,  1 ^2  ~ ^^0^1 


It  can  be  seen  from  this  equation  that  for  small  angular  misalignments 
this  particular  coefficient  will  experience  only  small  changes.  The  actual 
forms  of  the  transformation  equations  for  any  particular  test  item  will  depend 
on  the  magnitude  of  the  misalignment,  so  simple  algebraic  forms  may  result 
in  some  cases. 
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7,  Transformation  Sensitivity  Equations 

In  order  to  evaluate  the  magnitude  of  axis-offset  effects  on 
the  various  parameters  of  interest,  sensitivity  relations  can  be 
derived.  These  terms  can  then  be  evaluated  for  any  specific  case  to 
determine  whether  the  changes  in  parameters  are  significant  or  not. 

By  performing  this  calculation  prior  to  simulation  it  may  be  possible 
to  reduce  the  computing  time  for  some  configurations. 

Typical  results  are  shown  below  for  an  Intertia  component  and 
an  aerodynamic  coefficient.  A table  of  equations  for  other  parameters 
is  presented  in  Appendix  A. 


The  sensitivity  equation  for  the  inertia  component,  Ixx2> 
given  by  (for  a typical  angular  deviation) 


TT 


) Sin  2 J cos^e  - 2 cos  zS cos^e. 
- sinze  ( si’n cT  - 2^  cos  ) 


(IV-8a) 


The  corresponding  equation  for  a translational  deviation  is  given  by 

" ro  [zf  (l-^iosVcos^e)  - coi^e) 

+ h ( cos<rsrn  2^)]  (iv-8b) 


A similar  set  of  sensitivity  parameters  for  the  aerodynamic 
coefficient,  , is  given  by 

' " A,[-5i’n<5  (cos  € cosff  + z Siio2€  sin2(r 

V sih  6cos<r  ) f cost^cose  siio  2(r 

■=■  o 


2^, 


(IV- 9a) 

(IV- 9b) 


The  availability  of  sensitivity  equations  such  as  equations  (8) 
and  (9)  makes  possible  a determination  of  the  relative  magnitudes  of 
the  deviations  prior  to  performing  simulations. 
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SECTION  V 


RESEARCH  PLAN 


1.  Introduction 

The  preceding  sections  of  this  report  have  described  specific 
tasks  which  were  performed  to  provide  background  information  and  to 
generate  operational  tools  which  would  be  used  to  perform  more  detailed 
study  and  development  of  aerodynamic  models  to  be  used  in  the  analysis 
of  aeroballistic  test  data.  The  intent  of  such  a study  would  be  to 
improve  the  analysis  capabilities  for  conventional  munitions  but  also 
to  extend  the  envelope  of  applicability  of  ballistic  range  tests. 

In  this  section  a brief  discussion  of  a research  plan  directed 
toward  such  analysis  improvements  will  be  given.  Block  diagrams 
identifying  each  phase  of  the  project  will  be  discussed,  and,  to  the 
degree  possible  prior  to  actual  initiation  of  the  research  effort, 
the  methodology  to  be  employed  will  be  identified. 

2 . Simulation  Methodology 

A significant  part  of  the  aerodynamic  model  study  will  be  per- 
formed with  the  aid  of  a six  degree-of-freedom  computer  program  which 
solves  the  equations  of  motion  for  a complex  dynamic  system.  A basic 
program  written  at  Litton  Industries  (Reference  21),  the  M0D6DF  program, 
is  used  for  this  purpose.  The  program  is  structured  in  modular  form 
so  that  the  integration  algorithm  does  not  have  to  be  altered  for 
various  applications  of  the  program.  In  fact  alterations  of  the 
aerodynamic  models  contained  in  the  equations  of  motion  can  be  readily 
changed  without  affecting  the  remainder  of  the  program. 

A set  of  operational  subroutines  has  been  prepared  for  use 
within  the  basic  M0D6DF  structure.  These  subroutines  contain  the 
equations  of  motion  and  the  current  aerodynamic  models  which  are  being 
used  in  ballistic  data  analysis  at  Eglin  AFB.  The  equations  as  they 
appear  in  the  program  listing  are  shown  in  Appendix  B.  These  equations 
were  developed  by  General  Electric  and  are  shown  in  detail  in  Reference 
13.  The  current  aerodynamic  models  are  in  polynomial  form.  However, 
the  program  structure  is  such  that  any  analytical  model  could  be  used. 

The  standard  output  package  of  the  M0D6DF  program  has  been  modified 
to  be  corapatiblle  with  data  acquired  in  the  ballistic  range.  Output 
can  be  commanded  to  occur  at  the  downrange  positions  which  coincide  with 
data  stations  in  the  range.  This  simplifies  the  procedures  for  comparing 
simulated  and  free  flight  data.  Of  course,  time  increment  outputs  are 
also  available. 

An  additional  subroutine  containing  the  transformation  equations 
shown  in  Section  IV  of  this  report  has  also  been  Included  in  the  program. 
This  subroutine,  listed  in  /npendix  B,  permits  direct  accounting  for 
mass  offsets  in  the  simulation. 
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The  parameter  estimation  programs  at  Eglin  AFB  will  be  used  in 
the  data  analysis  phase  of  the  study.  These  programs  are  described 
in  Reference  13.  The  maximum  likelihood  method  of  estimation  is  used 
in  the  programs.  Thus,  the  aerodynamic  models  and  the  sensitivity 
parameters  which  involve  the  models  analytically  are  "hard  wired"  into 
the  program.  Spontaneous  changes  in  the  models  incorporated  into  the 
extraction  programs  are  not  possible.  This  implies  that  the  principal 
use  of  these  programs  will  be  in  performing  studies  with  the  existing 
models  to  determine  extraction-order  sensitivites  and  model  bias.  If 
significant  model  forms  which  have  been  proven  by  simulation  studies  do 
become  available,  the  extraction  programs  will  be  altered  to  accommodate 
them. 


Several  methods  of  parameter  estimation  were  discussed  in  Section 
II  of  this  report.  The  model  structure  determination  (MSD)  method 
appeared  to  have  considerable  promise  with  regard  to  the  present 
research  objectives.  However,  no  work  is  planned  in  this  area.  Rather, 
the  effort  will  be  directed  toward  improvement  of  the  existing  MLW  program. 

4.  Interactive  Methodology 

It  is  intended  that  most  of  the  analytical  studies  be  performed 
in  the  interactive  mode  at  the  graphics  terminal.  The  Tektronix  4014 
unit  is  a large  screen  terminal  with  extensive  graphics  capability.  A 
program  such  as  the  M0D6DF  simulation  program  can  be  run  from  the 
terminal  with  graphics  output.  A software  package  has  been  prepared 
which  will  permit  a frame  of  data  to  be  stored  and  then  recalled  in  any 
scale  and  to  any  part  of  the  screen.  This  software  allows  frames  to  be 
multiply  overlaid  so  that  simulation  runs  with  different  model  forms  can 
be  instantly  compared  with  each  other  and  with  test  data.  This  feature 
will  permit  rapid  identification  of  specific  characteristics  of  the 
various  models.  Of  course,  the  numerical  information  which  is  displayed 
graphically  on  the  screen  is  also  readily  available  for  detailed  study, 

5.  Synthesis  Methodology 

A block  diagram  of  the  interactive  computational  system  is  shown 
in  Figure  . Note  that  the  parameter  estimation  programs  are  actually 
installed  and  operational  on  the  Eglin  AFB  CDC  6600  computer  system. 

These  programs  will  be  run  by  remote  telephone  connection  on  the  graphics 
terminal  at  the  research  site.  This  procedure  eliminates  costly  instal- 
lation and  running  expenses  at  the  research  site. 

All  of  the  simulation  work  is  done  on  site  at  the  graphics 
terminal.  Since  this  operation  is  performed  within  Department  facilities, 
it  is  not  time  restricted. 

6.  Scope  of  Research 

The  overall  purpose  of  the  continuing  research  program  is  to 
provide  improved  aerodynamic  model  forms  for  use  in  ballistic  range  data 
analysis.  Three  major  tasks  can  be  identified  in  this  regard.  These  are: 
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(1)  Determination  of  the  effects  of  coefficient  extraction  order  on 
the  accuracy  of  results,  (2)  Aerodynamic  model  improvements  through  the 
use  of  alternative  model  forms,  and  (3)  Correlation  of  model  form 
selection  with  fundamental  aerodynamic  flow  field  phenomena.  Of  course, 
these  tasks  are  interdependent  and  cannot  be  separated  on  a phenomeno- 
logical basis,  but  the  methodology  to  be  employed  in  each  task  can  and 
will  be  discussed  separately.  It  should  also  be  pointed  out  that  the 
effort  will  be  directed  toward  studies  of  tail-finned  configurations  only. 

Apparently,  the  quality  of  the  parameter  estimations  performed  by 
the  AFATL  maximum  likelihood  program  is  affected  by  the  order  in  which 
extractions  are  made.  This  is  probably  related  directly  to  the  analytical 
model  which  is  chosen  to  represent  the  physical  system  in  the  estimation 
procedure.  In  the  present  research  no  effort  will  be  made  to  derive 
uniqueness  theorems  (or  lack  thereof)  for  the  estimation  procedure. 

Rather  the  method  will  be  applied  repetitively  for  various  models  and 
configurations  in  an  attempt  to  establish  model-related  trends.  In  this 
study  the  Aeroballistic  Research  Facility  Data  Analysis  System  (ARFDAS) 
developed  by  General  Electric  Company  for  AFATL  will  be  used.  The 
existing  aerodynamic  models  will  be  used  in  this  study  initially  since 
a considerable  supply  of  data  is  available  for  comparison  with  the 
present  analysis.  The  technique  to  be  employed  is  straight-forward; 
perform  the  parameter  estimation  procedure  systematically  on  different 
sets  of  coefficients  in  the  model,  compare  the  results  of  these  studies 
with  each  other  and  with  existing  data,  and  formulate  trends  in  the 
results.  Of  course,  systematic  manipulation  in  such  a manner  for  all 
combinations  of  coefficients  would  be  prohibitively  expensive,  so 
coefficients  which  have  historically  been  sensitive  to  such  manipula- 
tions will  be  isolated  for  the  study.  Certain  axial  force  coefficients 
and  all  rolling  motion  coefficients  are  prime  candidates  for  this 
treatment. 

The  benefits  to  be  gained  from  this  type  of  analysis  are  potentially 
great,  but  there  are  several  uncertainties  in  the  analysis.  Principal 
among  these  is  the  degree  to  which  the  effect  of  extraction  order  is 
dependent  on  the  model  form.  The  studies  will  be  done  only  with  existing 
models,  so  results  achieved  may  be  inconclusive  in  the  general  sense. 

The  statistical  nature  of  the  estimation  procedure  will  make  it  very 
difficult  to  determine  why  a particular  trend  occurs.  In  addition,  it 
may  turn  out  that  order  changes  may  produce  opposite  effects  in  different 
parameters.  If  such  occurrences  cannot  be  resolved  through  correlation 
with  the  actual  aerodynamics,  little  of  value  will  come  out  of  the  study. 
To  date,  no  comprehensive  investigation  has  been  performed  in  this  area, 
so  the  direct  method  of  analysis  is  a logical  choice.  In  any  case,  the 
expected  result  of  the  study  is  an  indication  of  the  type  and  degree 
of  interaction  which  exists  between  the  order  and  the  accuracy  of  the 
estimated  coefficients. 

The  remaining  two  tasks  associated  with  the  project  are  very 
closely  related.  Certainly,  the  aerodynamic  flow  field  produced  around 
the  flying  finned  body  directly  determines  the  aerodynamic  forces  and 
moments  on  the  body  which  dictate  its  motion.  The  analytical  model  of 
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these  aerodynamic  forces  and  moments  which  is  used  in  the  data  analysis 
algorithm  should  also  be  influenced  by  the  physics  of  the  flow  field. 
Reference  13  contains  the  analytical  models  which  are  currently  in  use 
at  AFATL.  For  many  cases,  these  models  yield  splendid  results.  However, 
there  is  still  some  uncertainty  in  the  models  for  some  cases.  For 
example,  in  one  recent  test  the  motion  of  the  test  item  was  very  well 
predicted  except  for  the  roll  orientation.  In  another,  the  axial  force 
estimates  were  apparently  in  error.  Thus,  some  model  improvements  are 
necessary,  but  more  importantly,  a systematic  method  for  selecting  a 
model  form  is  sought. 

Again,  a direct  method  of  analysis  is  suggested.  The  utility  of 
the  method  will  be  enhanced  by  the  use  of  computer  graphics  techniques 
which  permit  rapid  comparisons  of  a large  number  of  model  forms  and 
flow  field  effects.  A typical  model  analysis  procedure  can  be  outlined 
as  follows: 

(1)  Experimental  data  from  a test  of  a finned  body  are  stored. 

(2)  The  results  of  a parameter  estimation  run  for  that  case  using 
the  existing  aerodynamic  models  are  determined  and  stored. 

(3)  A 6D0F  simulation  of  the  test  item  motion  using  the  extracted 
coefficients  is  performed  and  stored  to  be  used  as  an  analytical 
base  1 ine . 

(4)  Simulations  using  variations  of  the  baseline  analytical  model 
are  performed,  and  the  graphical  results  are  stored  by  frame, 

(5)  The  extracted  simulation,  the  experimental  data  and  the  results 
using  new  aerodynamic  models  are  then  manipulated  interactively 
at  the  graphics  terminal  to  evaluate  the  influence  of  the  model 
changes  on  the  simulated  results.  Frequencies,  amplitudes  and 
phase  information  are  compared  for  the  various  cases  to  determine 
model  change  effects. 

The  present  research  in  model  form  will  not  be  confined  to  seeking 
higher  order  effects  in  the  conventional  expansion  of  aerodynamic  forces 
and  moments.  Some  forms  of  aerodynamic  data  may  suggest  unexpected 
types  of  non-linearities  which  cannot  be  treated  in  the  usual  manner. 

An  example  of  this  type  which  has  already  been  observed  is  a drag  non- 
linearity which  is  not  symmetrical  and,  thus,  cannot  be  modeled  with 
an  even-power  expansion.  Other  new  forms  may  also  be  dictated  by  the 
unique  aerodynamic  regimes  in  which  some  of  the  new  munitions  operate. 
High  angle  of  incidence  is  common  in  these  casses,  and  the  large 
separated  regions  and  strong  vortices  produce  unusual  behavior.  Because 
the  simulation  programs  to  be  used  in  the  study  are  modular,  the  new 
aerodynamic  model  forms  which  may  be  dictated  by  these  flow  fields  can  be 
easily  incorporated  Into  the  programs.  Even  for  radical  model  forms, 
the  procedure  to  be  followed  will  be  as  outlined  above.  It  should  be 
understood,  however,  that  these  new  models  will  not  be  incorporated  into 
the  parameter  estimation  programs  until  the  selected  forms  have  been 
shown  to  be  superior  for  several  different  configurations. 

Correlation  of  flow  field  information  with  aerodynamic  model  form 
will  be  accomplished  by  analysis  of  wind  tunnel  data,  flow  field  pictures 
and  ballistic  test  data.  The  emphasis  here  will  be  on  attempting  to 
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detect  specific  features  of  the  flow  field  which  manifest  themselves 
in  certain  aspects  of  the  model  form.  As  an  example,  an  induced  rolling 
moment  may  be  attributable  to  body-fin  interference  which  could  be 
described  in  the  aerodynamic  model  in  terms  of  certain  aerodynamic  or 
geometric  parameters.  The  basic  procedure  here  will  be  to  seek  trends 
in  behavior  and  to  relate  those  trends  to  the  model  form. 

The  proposed  method  for  dealing  with  the  aerodynamic  model  forms 
is  direct.  Such  an  approach  may  be  expected  to  discover  major  deviations 
from  conventional  model  forms.  It  may  not,  however,  provide  detailed 
information  on  subtle  flow  field  effects  on  the  models.  In  such  cases 
where  obvious  discrepancies  occur,  additional  perturbations  of  the 
evaluation  procedure  may  be  necessary.  Although  the  nature  of  these 
possible  deviations  cannot  be  predicted  now,  it  is  felt  that  the  basic 
tools  of  analysis  are  sufficiently  flexible  to  allow  for  them  in  some 
manner.  Systematic  variation  of  model  form  within  the  simulation 
programs  would  be  possible,  and  the  effects  of  these  variations  on 
specific  parameters  could  be  investigated  using  the  frame  overlay 
capability  of  the  graphics  terminal. 


19 


SECTION  VI 


CONCLUSIONS 


In  the  foregoing  report  a research  plan  has  been  described, 
and  the  basic  computational  tools  required  to  implement  the  plan 
have  been  developed.  The  study  has  indicated  that  model  improvements 
are  necessary  if  the  envelope  of  applicability  of  the  ballistic  range 
as  a test  facility  is  to  be  expanded  to  include  airplane- 1 ike  test 
items.  The  approach  to  improving  these  models  has  been  outlined  here. 
The  present  study  has  been  concerned  with  exploring  the  possibilities 
of  research  in  this  area.  An  effort  to  perform  model  development 
has  not  been  made. 
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FIGURE  2.  COMPUTATIONAL  FACILITIES 


TRANSFORMATION  EQUATIONS 
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